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Abstract: 

Absolutely  calibrated,  time-resolved  spectral  intensity  measurements  of  soft  x-ray  emission  from 
laser-irradiated  polystyrene  targets  are  compared  to  radiation-hydrodynamic  simulations  that 
include  our  new  postprocessor,  Virtual  Spectra.  This  new  capability  allows  a  unified,  detailed 
treatment  of  atomic  physics  and  radiative  transfer  in  non-LTE  conditions  for  simple  spectra  from 
low-Z  materials  as  well  as  complex  spectra  from  high-Z  materials.  The  excellent  agreement 
(within  a  factor  ~  1 .5)  demonstrates  the  powerful  predictive  capability  of  the  codes  for  the 
complex  conditions  in  the  ablating  plasma.  The  absolute  spectral  intensity  measurements  were 
made  in  the  XUV  region  (hv  ~  0.1  -1.0  keV)  with  transmission  grating  spectrometers  with  good 
time  resolution  (t  ~  0.3  ns)  and  moderate  spectral  resolution  (E/8E  ~  10).  Comparison  to  data 
with  high  spectral  resolution  (E/8E  ~  1 000)  emphasizes  the  importance  of  including  radiation 
coupling  in  the  quantitative  simulation  of  emission  spectra. 

The  knowledge  of  absolute  intensity  of  XUV  and  x-ray  radiation  created  by  the 
laser  irradiation  of  solid  targets  is  crucial  for  issues  in  inertial  confinement  fusion  (IGF) 
ranging  from  illumination  sources  for  backlighter  diagnostics1  to  efficient  sources  to 
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drive  implosions.2  The  correct  estimation  of  the  role  of  radiation  transfer  in  the  evolution 
of  the  plasma  depends  in  a  large  measure  on  the  intensity  of  the  emitted  radiation.  This  is 
especially  true  for  indirect-drive  ICF  which  utilizes  an  envelope,  or  hohlraum,  around  the 
fuel  pellet  to  convert  laser  light  into  smoothed  XUV  radiation  which  ablates  and 
implodes  the  pellet.  Typically,  hohlraums  are  made  from  Au  and,  accordingly,  the 
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absolute  XUV  emission  from  laser-irradiated  Au  has  been  studied  extensively. 

However,  most  pellet  designs,  whether  for  indirect-drive  or  direct-drive  ICF,  incorporate 
a  polystyrene  (CH)  shell  to  contain  the  DT  fuel.  This  central  role  of  CH  means  that  a 
comparison  of  absolute  measurements  of  its  intensity  from  a  laser  target  with  simulations 
of  absolute  emissivity  is  vital.  Numerous  spectroscopic  studies  of  laser- irradiated  low-Z 
targets  such  as  CH  have  been  published4  but,  to  the  best  of  our  knowledge,  an 
experimentally  benchmarked  prediction  of  absolutely-calibrated,  spectrally-resolved 
emission  from  CH,  relevant  to  ICF  purposes,  has  not  been  published.  In  this  paper,  we 
compare  the  emission  from  planar  CH  targets  measured  with  absolutely  calibrated 
transmission  grating  spectrometers  to  simulations  performed  with  a  new  post-processor 
(Virtual  Spectro,  or  VS).  We  report  agreement  between  the  predicted  and  measured 
values  for  the  absolute  intensity  better  than  a  factor  of  ~1.5  in  the  soft  x-ray  region  (hv  ~ 
0. 1  - 1 .0  keV).  This  agreement  provides  confidence  in  the  radiation  hydrodynamic 
simulations  used  at  U.S.  Naval  Research  Laboratory  (NRL)  for  ICF  pellet  design,  in 
which  the  opacities  are  calculated  by  the  STA  model5  and  the  non-LTE  effects  are  treated 
with  the  RADIOM6  model.  The  post  processor  VS  goes  beyond  STA  and  LTE 
approximations,  it  is  a  detailed  fine- structure  radiative  collisional  model  with  radiation 
coupling. 


2 


Absolute  comparisons  of  spectral  intensities  are  seldom  made,  because  they  are 
difficult.  On  the  experimental  side,  everything  has  to  be  calibrated:  the  laser  irradiance, 
on  the  targets,  all  the  components  of  the  spectrometers  -  slits,  fdters,  gratings,  detectors, 
etc...-  as  well  as  geometrical  factors,  like  solid  angles  and  distances.  On  the  theoretical 
side,  the  spatial  and  temporal  histories  of  temperature  and  density  of  the  plasma  need  to 
be  accurately  reproduced,  reliable  atomic  physics  models  are  necessary  for  the  population 
of  the  emitting  energy  levels  and  for  transition  probabilities,  and  radiative  transfer  and 
coupling  to  the  various  atomic  species  has  to  be  taken  into  account. 

Radiation-hydrodynamic  simulations  of  laser-irradiated  targets  must  include  a 
large  number  of  phenomena  to  provide  reliable  descriptions  of  experimental  conditions. 
Work  at  NRL  has  focused  on  developing  a  code  suite  to  analyze  direct  drive  schemes  for 
ICF.  The  hydrodynamic  code,  FAST,7  is  a  multidimensional  code  which  has  been 
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parallelized  to  enable  faster  operation  for  2D  or  3D  runs.  It  uses  flux  corrected  transport 
to  provide  a  highly  sophisticated  adaptive  gridding  algorithm  that  efficiently  captures  the 
sharp  gradients  inherent  in  these  physical  systems.  FAST  works  with  the  quasineutral 
approximation  to  treat  the  plasma  as  a  single  fluid,  although  distinct  electron  and  ion 
temperatures  can  be  maintained  in  the  runs.  Laser  absorption  is  described  by  inverse 
bremstrahlung  with  induce  spatial  incoherence  smoothing  (ISI).  Electron  thermal 
conduction  follows  classical  flux  limited  Spitzer-Hann  model  with  quantum  corrected 
conductivity.9  The  various  equations  of  state  for  the  target  materials  are  tabulated,  and  a 
multigroup  radiation  diffusion  model  with  variable  Eddington  factor  is  used  for  radiation 

transport^.  The  necessary  opacities  are  interpolated  for  each  cell  from  tables  calculated 
off-line  for  a  dense  grid  of  values  for  temperature,  density,  and  photon  energy.  In  the 
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present  case,  these  tables  were  obtained  with  the  STA  model  in  local  thermodynamic 
equilibrium  (LTE),  and  corrected  for  non-LTE  effects  with  the  RADIOM  model.  FAST 
can  also  use  tables  from  SCROLL,10  a  non-LTE  version  of  STA,  but  this  was  not 
necessary  for  this  paper,  due  to  the  relatively  low  emissivity  of  CH.  The  FAST  code  has 
been  successfully  tested  against  primarily  hydrodynamic  effects  (e.g.  Rayleigh-Taylor 
instability  growth11,  Richmeyer-Meshkov  instabilities12).  Now  detailed  predictions  for 
radiation  effects  have  required  development  of  a  new  post-processor. 

Like  FLY14  or  TRANSPEC, 15  Virtual  Spectro  (VS)  16  is  a  post-processor  which 
uses  temperatures,  densities,  average  charge  state,  and  velocities  as  a  function  of  time  and 
position  as  computed  by  a  ID  or  2D  hydrodynamic  code.  The  hydrodynamic  code  dumps 
from  FAST  can  be  given  in  either  Lagrangian  mode,  Eulerian,  or  arbitrary  Lagrangian- 
Eulerian  (ALE)  mode.  VS  can  process  a  pure  material  or  each  component  of  a 
compound.  First,  using  the  mass  density  and  electron  temperature  for  each  cell,  VS 
computes  the  relevant  specific  electron  density  for  the  simulation.  The  atomic  structure 
data  (taken  from  published  values  or  computed  with  the  HULL  AC  10  code  suite)  is  then 
used  to  compute  time  dependent  atomic  populations,  outgoing  flux  (Jv),  and  opacity  (Kv) 
with  a  second  order  advection  scheme  in  the  non  ALE  mode,  though  it  has  been  checked 
that  the  quasi-steady  state  approximation  is  valid  for  the  present  study.  Radiation 
transport  in  resonance  lines,  taking  into  account  line  broadening,  is  coupled  to  kinetic 
equations  using  a  mix  of  the  accelerated  lambda  iteration  scheme  and  the  partial 
linearization  method,  see  Fig.  1.  The  final  outputs  are  used  to  create  synthetic  spectra 
with  temporal,  spatial,  and  energy  resolution  where  opacity  effects  are  accurately 
accounted  for.  Experimental  setup  (including  integration  along  the  line  of  sight  at  a 
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particular  view  angle,  instrumental  spectral  resolution,  and  time  integration)  can  be 
incorporated  in  the  final  output  spectra. 

Experiments  to  test  the  predictions  of  FAST  and  VS  were  carried  out  at  the  Nike 
laser  facility.  The  Nike  laser17  is  a  high  power  (-10 14  Wcm'2)  KrF  system.  Nike  operates 
at  248  nm,  has  a  relatively  large  bandwidth  (1-2  THz),  and  provides  very  uniform  beam 
intensity  profiles  (5I/I  <  0.2  %)  on  a  500pm  focal  spot  due  to  the  use  of  the  ISI  smoothing 
technique.  For  the  results  in  this  paper,  standard  40  pm  thick  planar  CH  targets  were 
placed  at  target  chamber  center  and  irradiated  with  a  3.5  ns  pulse  at  an  incident  power 
density  of  1.13x10  ~  Wcm"  .  The  prepulse  intensity  was  less  than  a  factor  10  below  the 
level  of  the  main  pulse.  The  normal  to  the  target  surface  was  aligned  along  the  average 
laser  propagation  direction  for  the  laser  beams. 

Soft  x-rays  emitted  from  the  target  plasma  were  recorded  with  two  transmission 
grating  spectrometers  (TGS)  and  a  grazing  incidence  spectrometer  (GIS).  These  TGS 
were  installed  to  view  the  target  at  two  widely  separated  viewing  angles  to  determine  the 
angular  dependence  of  the  soft  x-ray  emission.  The  ‘normal’  viewing  spectrometer  had  a 
line  of  sight  at  an  angle  of  14°  to  the  target  normal  and  the  ‘wide’  viewing  spectrometer 
had  a  line  of  sight  at  42°  to  target  nonnal.  The  TGS  provided  time-integrated  and  time- 
resolved  spectra  (St  -  0.3  ns)  with  moderate  spectral  resolution  (87JX  -10-20)  and 
accurate  determination  (ai/I  <  15%)  of  the  intensity  (W  stcr  1 ).  The  main  optical 
elements  of  both  TGS  were  a  2500  1/mm  gold  transmission  grating,  an  adjustable 
entrance  slit,  a  free-standing  1800  A  thick  A1  foil  transmission  filter,  a  Si  photodiode 
array,  and  a  phosphor-coupled  CCD  camera.  The  sensitivity  calibrations  of  the 
transmission  filters,  transmission  gratings,  and  photodiode  arrays  were  performed  at  the 
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X24C  beam  line  of  the  National  Synchrotron  Light  Source  at  Brookhaven  National 
Laboratories.19 

In  this  paper,  we  show  only  results  from  the  ‘normal’  viewing  spectrometer  due 
to  its  higher  spectral  resolution  and  because  its  line  of  sight  is  more  closely  aligned  to  the 
target  normal.  Observations  that  included  the  wider  viewing  spectrometer  verified  that 
the  variation  of  absolute  intensity  with  view  angle  is  not  significant  for  the  slight  off- 
normal  view  angles  in  the  case  of  CH  targets. 

The  GIS  provided  high  spectral  resolution  (8AA  -1000)  but  was  time-integrated 
and  not  absolutely  calibrated.  The  GIS  included  a  cylindrically  curved  focusing  mirror,  a 
40  pm  wide  entrance  slit,  a  1200  1/mm  gold  coated  concave  grating  with  a  2°  blaze  angle 
and  a  3  m  radius  of  curvature.20  Time-integrated  spectra  were  recorded  with  glass 
photographic  plates.  We  calibrated  the  intensity  response  of  the  film  plates  from  the  ratio 
of  the  line  intensities  using  the  standard  procedure  of  multiple  shot  exposures.  The 
grazing  incidence  spectrometer  viewed  the  target  at  an  angle  of  42°.  The  spatial 
resolution  across  the  target  at  the  photographic  plate  was  -50  pm/mm. 

We  have  post-processed  hydrodynamic  simulations  computed  in  the  one 
dimensional  mode  of  FAST  (FAST  ID)  using  400  spatial  cells  and  98  non-equidistant 
energy  groups  to  provide  adequate  resolution  of  the  carbon  lines  while  allowing  tractable 
computation  times.  The  laser  pulse  time  history  and  peak  intensity  (1.13x10  ~  WcnT  ) 
were  taken  from  the  laser  diagnostics  for  the  observed  shots..  The  FAST  output  was 
collected  at  80  dump  times.  In  VS,  the  atomic  structure  was  taken  from  the  literature  for 
H-like  and  He-like  resonant  lines  and  for  satellites  of  a  and  p  lines.  In  all  ions, 
hydrogenic  levels  have  been  added  to  complement  the  database  up  to  n=15.  Stark 
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broadened  line  profiles  using  tabulated  microfield  distribution  were  used  for  the  resonant 
lines  up  to  n=4.22  As  mentioned  previously,  only  resonant  lines  were  considered  here 
due  to  the  comparatively  low  intensity  for  the  satellite  lines.  For  the  comparison  with 
experiment,  the  synthetic  spectra  are  convolved  with  (gaussian)  instrumental  functions  of 
widths  appropriate  for  each  instrument. 

Normal  incidence  synthetic  spectra  in  the  range  250-600  eV  and  time-  or  spectral- 
integrated  values  are  compared  to  the  absolute  intensity  measurement.  Figure  2 
demonstrates  that  the  simulation  reproduces  correctly  the  rise  of  the  emission.  Note  that 
the  simulations  in  both  Fig.  2  and  3  are  multiplied  by  a  factor  1.5  to  emphasize  the 
agreement  on  the  shapes.  After  the  laser  pulse  ends,  the  simulation  has  a  slower  decay  of 
the  emission  compared  to  the  measurements.  This  may  be  due  to  2D  cooling  effects  not 
included  in  the  ID  simulation.  Figure  3  compares  the  observed  values  for  the  absolute 
time-integrated  spectral  intensity  (J/ster/A)  to  the  final  output  from  VS.  The  two  absolute 
spectra  agree  very  well.  There  is  an  apparent  wavelength  shift  of  the  satellite  lines  at  38 
A  .  Several  lines  are  blended  there  due  to  the  poor  resolution.  A  slight  error  in  the 
temperature  during  the  plasma  expansion  after  the  laser  pulse  could  change  the  intensity 
ratios  and  cause  this  apparent  shift.  The  -25%  discrepancy  in  the  intensity  ratio  of  He- 
like  to  H-like  emission  may  also  be  due  to  such  an  error. 

Figures  4  and  5  compare  a  high-resolution  time  integrated  spectrum  recorded  on 
film  by  the  GIS  with  corresponding  VS  results.  They  also  show  simulations,  with  and 
without  coupling  to  the  radiation  field,  of  the  level  populations  up  to  n=4  for  hydrogen¬ 
like  and  helium-like  carbon.  All  spectra  in  both  figures  have  been  normalized  to  the  peak 
of  the  Lyp  line.  It  can  be  seen  that  the  inclusion  of  line  coupling  for  computation  of 


7 


populations  causes  a  dramatic  change  in  the  simulated  spectra.  The  unobserved  features 
around  27  A  of  Fig.  4  disappeared,  and  the  ones  at  34  A  and  42  A  are  much  smaller, 
although  still  present.  Also  the  level  of  the  continuum  below  26  A  agrees  better.  Most  of 
the  remaining  discrepancies  are  in  the  line  intensity  ratios,  which  depend  on  the  relative 
abundances  of  the  H-like  and  He-like  ions.  These  differences  can  be  attributed  to  slight 
errors  in  the  simulated  temperatures  from  FAST,  or  to  transitions  in  the  free-bound 
continuum  not  coupled  through  photo-ionization,  or  to  the  limitations  of  the  atomic 
model.  It  is  clear,  however,  that  radiation  transport  and  coupling  in  the  rate  equations  of 
atomic  levels  play  an  important  role. 

All  in  all,  the  agreement  between  simulation  and  experiment  shown  in  Fig.  2  and 
3  is  remarkable,  in  view  of  the  difficulty  of  the  calibration,  and  the  complexity  of  the 
simulation.  This  has  important  implications:  because  the  underlying  data  for  Virtual 
Spectro  is  the  output  from  the  FAST  ID  code,  the  success  of  the  final  product  implies  that 
the  radiation  hydrodynamics  code  suite  is  predicting  the  temperature  and  density  of  the 
emission  region  accurately.  The  radiation  present  in  the  physical  system  does  not 
dominate  the  evolution  of  the  system.  The  STA  tables  and  RADIOM  models  provide 
fairly  good  accuracy,  although  not  the  detailed  spectral  infonnation  created  by 
Virtual  Spectro.  These  new  results  represent  a  significant  advance  in  the  ability  to 
accurately  describe  conditions  in  laser-irradiated  targets.  We  plan  to  address  the 
comparison  of  spectra  of  high-Z  elements,  which  is  much  more  difficult,  in  the  near 
future. 

This  work  was  supported  by  the  U.  S.  Department  of  Energy. 
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Figure  Captions 

1 .  Radiation  trapping  -  between  different  zones  -  increases  excited  states 
populations  (thus  local  emissivity),  but  decreases  the  outgoing  x-ray  flux. 

2.  Spectrally-integrated,  time  resolved  data.  Comparison  of  outputs  from  Virtual 
Spectro  (dashed  line)  and  from  absolutely  calibrated  transmission  grating 
spectrometer  (solid  line  and  closed  circle).  The  simulated  spectral  intensity  was 
multiplied  by  1.5. 

3.  Time  integrated,  spectrally  resolved  data.  Comparison  of  outputs  from  Virtual 
Spectro  (dashed  line)  and  from  absolutely  calibrated  transmission  grating 
spectrometer  (solid  line  and  closed  circle).  The  simulated  spectral  intensity  was 
multiplied  by  1.5. 

4.  Time  integrated  spectrum.  The  comparison  of  simulation  without  coupling 
between  resonant  lines  (dashed  line)  to  high  resolution  spectrum  from  grazing 
incidence  spectrometer  (solid  line)  shows  an  overestimate  of  the  intensities  of 
many  spectral  lines,  particularly  satellite  lines.  Each  spectrum  has  been 
nonnalized  to  the  Ly[!  intensity  of  hydrogen-like  carbon. 

5.  Time  integrated  spectrum.  The  comparison  of  simulation  with  coupling  between 
resonant  lines  (dashed  line)  to  high  resolution  spectrum  from  grazing  incidence 
spectrometer  (solid  line)  shows  an  improved  estimation  of  the  intensities.  Each 
spectrum  has  been  normalized  to  the  Ly[!  intensity  of  hydrogen-like  carbon. 
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Figure  2 
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